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CHAPTER  1 
INTRODUCTION 


1.  Introduction  and  Summary 

O/w- 

*~^S>  In  this  dissertation  w e-study  stochastic  piocesses^hich  arise  out 

jjb, 

of  a single  Gl/G/l  queueing  system.  W.e  are>concern«i  herewith  the 
weak  convergence  of  each  process  conditioned  on  an  event  whose  probability 
is  converging  to  zero  as  time  becomes  large, 

\ 

The  basic  processes  we  are  going  to  discuss  are  {W^:  n £ 0}, 

{ ©( t) : t S 0}  and  {W(t):  t S 0),  where  represents  the  time  which 

customer  numb  * n must  wait  before  being  serviced,  o(t)  is  the  number 
of  customers  present  in  the  system  at  time  t,  and  w(t)  denotes  the 
workload  of  the  server  at  time  t. 

The  setup  for  our  problems  is  the  classical  one.  in  a single  Gl/G/l 

queueing  system,  customer  number  0 arrives  at  time  t^  = 0,  finds  a free 

server,  and  experiences  a service  time  v^.  The  n**1  customer  arrives  at 

time  t and  experiences  a service  time  v . Let  the  interarrival  times 
n n 

t -t  , k u , n 2 1.  We  assume  the  sequence  of  random  vectors 
n n-1  n 

{( v , u ) : n 2 1}  is  independent  and  identically  distributed  (i.i.d.). 
n-1  n 

Let  E(u^)  = X 1 and  Efv^}  = |i  \ where  0 < X,  p < ®.  In  addition,  we 

2 2 

shall  always  assume  that  E{u  } + E{v  } < ® and  that  the  deterministic 
J n n 

system  in  which  both  u and  v are  degenerate  is  excluded.  The  natural 

n n 

measure  of  congestion  for  this  system  is  the  traffic  intensity  p = X/p. 

In  our  study  we  shall  consider  systems  in  which  p is  greater  than,  equal 


to  or  less  than  unity. 


If  we  define  5 = v , - u , then  W can  be  defined  recursively 

n n-1  n n 

by  the  relations 


*0  " °’ 


W„+l  = [\  + W’  ”*°- 


Now  {§  : n S 1}  is  a sequence  of  i.i.d,  random  variables  with  E^}  = 

2 2 2 

and  E{5  — |i^t=a  > 0 < o“  < «,  Forn  the  random  walk  (S^:  n £ 0}  by  setting 

S.  = 0 and  S = £+•>.  •+|,n^l«  It  is  easy  to  show  by  induction  that 
0 n 1 n 


Wn  = maxfP^  - Sk  :k  = 0,  1,  . n),  r ;»  0 

Next,  let  N bo  the  hitting  time  of  the  set  (-“,  x]  by  the  random 
walk: 

N ~ inffn  > 0:  S ^ x}  , 
x 1 n 1 ! 

where  the  infimum  of  the  empty  set  is  taken  to  be  +“.  Set  N = Nq.  In 
the  queueing  context  of  is  the  number  of  customers  served  in  the  first 
busy  period. 

What  sort  of  limiting  results  would  v-e  hope  to  provide  for  {Wn:  n £ 0} 

and  the  other  processes?  One  would  like  to  have  results  of  convergence  in 

4 

distribution.  For  example,  does  there  exist  a random  variable  S , whose 
distribution  is  proper,  s"ch  that 

lim  PfS  i x I n < N < "]  = P{S^  £ xl 

ir>oo  1 n 1 

for  each  x £ 0 when  P>  1 (i.<-:.,  > 0)?  We  observe  that  Wn  = 

4 

on  {N  > n}.  Thus  S'  can  be  thought  of  as  a limit  law  for  the  waiting 
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time  process  in  the  first  busy  period  cf  a general  single  server  queue 
with  P > 1. 

Let  [jc]  be  the  greatest  integer  m x.  Define  the  random  function 

X Ly 
u 

xn(  t)  , ostsi, 

an5 

where  a is  a constant  to  be  specified  later.  Then,  we  would  like  to 
know  whether  one  car.  obtain  a functional  central  limit  theorem  as 

(X  | n < N < ®)  =»  X , 

where  X is  certain  random  function.  The  symbol  =*>  means  weak  converg- 
ence in  the  Skorohod  topology  in  D[0,lj  = D;  see  BILLINGSLEV  (1968), 
p.  137,  for  further  discussion  of  this  result. 

The  limiting  random  function  X that  we  shall  encounter  in  our  study 
can  be  identified  in  terms  of  Brownian  motion  B.  The  limit  is  either 
Brownian  meander  or  Brownian  excursion  depending  upon  whether  P = 1 or 

P / 1. 

Brownian  meander,  W+,  is  first  i eccgnized  in  BELKIN  (1972),  Theorem 
3.1.  Let  {B( t)  . t ^ 0}  be  the  standardized  Brownian  motion.  Let 
T1  * sup{t  e [0- 1 3s  B(  t)  = 0).  Set  = 1 - Then 


w+(t) 


|B(x1  + tA^I  / A^. 


c <;  t < i . 


It  is  s continuous,  non-homogen eons  Markov  process  and  has  transition 
density  given  by 
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(1.1)  P{W+(t)  t dy)  =r  P+(0, 0,  t, y)dy 

- A expp-jlNlf-^-jdy 

\ 2t/  '.(i_t)2/ 

for  0 < t £ 1 and  y > 0;  for  0 < s < t £ 1,  and  x,  y > 0 

(1.2)  P{W+(tJ  * dy  I W+(s)  = x) 

n|( 

= p+(s,x,  t,y)dy  = g(t-s,x,  y)  

|n|( 

2 2 

where  g(t,x,y)  = (2flt)  ^ Cexp(--^~“)-  exp(  ~)1  and 

ix  2 

|n|(x)  = (l/it)  j exp(-  -!~)du,  X > 0 . 

For  the  derivation  of  this  transition  density,  see  BELKIN  (1972),  p.  61. 
Note  that  for  x £ 0, 

x 

P{W+(l)  £ x)  = j R(y)dy  , 

0 

where  R(x)  = x exp(-  x2/2),  x ^ 0,  is  the  density  cf  the  Raleigh  distri- 
bution. 

Another  process  which  we  wi‘i  be  interested  in  is  the  Brownian  excur- 
sion, W+,  Let  i0  ~ inf{t  a 1:  w(t)  =0).  Set  A„  = t - i and 

»o(t)  = Ib^  +■  tA2 ) | / Ag  > H s t < 1. 

Brownian  excursion  is  also  a continuous,  non-homogeneous  Markov 
process  witn  transition  densitj  given  by 


-h- 


(1.4) 


P('wt(  t)  e dyj  = pn(0,0,t,y)dy 


o 2 

2y  e^p 


(-  W 

\ 2t(i-7T/ 


\&3(i-t)3 


for  0 < t < 1 and  y > 0;  for  0 < s < t < 1 and  x,  y > 0 


(1.5)  P(wJ(t)  e dy  | vj(s)  = x) 


= Pq(s,x,  t,y)dy  = g(t-s,  x,y)|^5fj 


3 y exp 

-s\2 


(-  irhr) 


x exy 


/ __x \ 

\ (2(1-8)/ 


dy  ; 


see  ITO-McKEAN  (1965),  p.  7 6,  for  this  result.  Clearly,  the  distribution 
of  W+(0)  and  Wq(i)  are  degenerate  at  0. 

In  Chapter  2 we  shall  show  that,  if  P £ 1,  the  fin ite-u.lraeis tonal 
distribution,*  (f.d.d.’s)  of  (X  | n < N < “)  converge  weakly  to  those  of 


W0* 


Because  of  the  degeneracy  at  point  t = 1,  the  limiting  behavior 
v W | n < N < “)  is  not  revealed.  When  P > J,  the  existence  of 

’Ct 

lim  (W  | n < N < ») 
n->®  r 


will  be  established  and  it,  distribution  function  will  be  exhibited. 
This  is  carried  out  in  Chapter  3» 

In  Chapter  4 we  restrict  our  attention  to  the  case  that  P = 1. 

We  establish  results  for  the  workload  process,  W(t),  and  the  queue- 
length  process,  y l).  Results  for  certain  multiple  channel  queues  are 
also  given. 
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{ 

j To  close  this  section,  we  make  a comment  regarding  the  notation  in 

i 

i this  thesis.  The  basic  notation  adhered  throughout  the  paper  is  devel- 

i 

oped  in  Chapter  1 and  the  first  sections  of  Chapters  2 and  4.  Bracket  , 

[ ],  refers  to  the  bibliography. 

2.  Related  Work 

The  history  of  conditioned  limit  results  for  various  stochastic  pro- 
cesses does  not  extend  over  the  past  twenty  years.  The  existence  of  con- 
tioned  limit  distributions  in  both  discrete  an'*  continuous  time  Markov 
chains  is  extensively  discussed  by  MANDEL  (1959,  I960),  OARROCH  and  SENETA 
(1965,  1967),  SENETA  and  VERi -TONES  (l967). 

Results  on  conditioned  random  walks  have  been  obtained  by  DffASS 
and  KARLIN  (1963),  LIGGETT  (1968),  DALEY  (1969),  BELKIN  (1970,  1972), 

KAIGH  (1974),  and  IGLEHART  (l974a,b). 

For  the  m/g/i  queueing  system  KYPRIANOU  (l97l)  first  established 
the  existence  of 

lim  PfW(t)  s x I t < T(a)<  «},  0 < x < °° 

t-»°=  ’ 

where  T(a)  = inf{t  : W(t)  = 0,  W(0)  = a)  for  a /=  0.  It  is  further 
. demonstrated  by  KENNEDY  (1974)  that  as  time  becomes  large  w(t),  suitably 

J 

scaled  and  normed,  converges  to  Brownian  excursion  or  Brownian  meander 
depending  upon  whether  P / 1 or  P = 1, 

The  first  conditioned  limit  result  for  the  Gl/G/i  queue  is  contained 
in  IGLEHART  (l974a  j.  In  [17]  he  established  the  following  result. 


1 


(1.6)  THEOREM.  [IGLEHART  (1974a),  Theorem  3.4].  _If  u = 0,  0 < a2  < «, 

E{  | ) < «o,  and  | nonlattice  or- integer-valued  with  span  1 , then 

(X  I N > n)  =$>  W+. 
n 

The  third  moment  condition  was  later  removed,  see  DURRETT  (1976). 

When  < 0,  which  corresponds  to  a stable  queue,  the  random  walk 

(Sn:  n S 0}  is  strongly  attracted  to  the  origin  and  the  limit  results  in 

2 

this  case  do  not  depend  only  on  and  o but  on  the  entire  distribu- 
tion of  Assuming  that  the  distribution  of  | satisfies  the  follow- 

ing coiiditions: 


(1.7) 

- £ p1  <0; 

(1.8) 

S(s)  = E exp(s  • g ) converges  for 
for  some  a > 0; 

real  se[0,a), 

(1.9) 

Q(s)  attains  its  infimum  at  a point 
where  &(t)  = Y < 1 and  9’(t)  =0; 

t,  0 < t < a, 

and 

(1.10) 

if  is  lattice,  then  P{g^  = 0)  > 

o, 

IGLEHART  (1974b)  established  the  following  result: 

(1.11)  THEOREM.  [IGLEHART  (1974b),  Theorem  2.3%  If  conditions  (.’.7)  - 
(1.10)  hold  and  u ^ 0,  then 

< “ y~n  ) 

<1*12>  E{exp(-uWn)|N  > n)  = expj  £ — [ E{exp(-uS*)  )-l]  j 


^ f(u)  . 


W*‘ VW?vs*«*B?aa*»iw«*»«fc< 


causer 


The  distribution  function  of  the  limit  random  variable  S*  is 


later  given  by  VERAVERBEKE  and  TEUGELS  (1975).  Let  U(^,  x)  =^2^  PfS^ 


x,  N > n}^  and 
(1.13) 


rx 

V(x)  = 1 - e ^ + T \ U(”,  x-y)e  TJ  dy. 

•>0+  { 


Then  under  the  same  conditions,  (1.7)  — (1.10),  it  is  shown  in  [303,  p.  283, 
that 


(1.14) 


lira  P{W  * x|n  > n) 
n-»oo  n 


= exp 


n=l 


y*n 

-i—  pf > 0} 
n n 1 


V(x), 


Acknowledgment:  We  close  this  introductory  chapter  by  a~  know-lodging  a 
heavy  debt  to  the  two  above-mentioned  papers  of  IGLEHART  (1974a,  b)  on 
conditioned  limit  theorems  for  random  walks # 
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CHAPTER  2 


LIMITING  DIFFUSION  FOR  RANDOM  WALKS  WITH  DRIFT 
CONDITIONED  TO  STAY  POSITIVE 


Introduction  and  Probability  Space 

Let  k S I]  be  a sequence  of  independent,  identically  distributed 

2 2 

(i.i.d.)  random  variables  with  Ef^)  = p^  and  E{|j  “ } = o , 

2 

0 < o < <o4  Form  the  random  walk  {S^:  n S 0}  by  setting  S^  = 0 and 

S =4,  +••*+£,  nSl.  Let  N be  the  hitting  time  of  the  set 
n *1  bn  x 

(-», x]  by  the  random  walk: 

N = inf {n  >0:  S £ x]  , 

where  the  infimum  of  the  empty  set  is  taken  to  be  +“.  Set  N = NQ. 

Throughout  this  chapter,  we  shall  assume  that  the  distribution  of 
4X  satisfies  the  following  conditions 

(2.1)  px  4 0; 

(2.2)  e(s)  = E exp(  £ ) converges  either  for  real  se[0,  a), 

if  Pj  < 0,  or  for  real  /4e(-a,0],  if  > 0,  for  some  a > 0; 

(2.3)  &(*)  attains  its  infimum  either  at  a point  t,  if  Pj  < 0, 
or  at  a point  -t,  if  p^  > 0,  where  0 < T < a,  and 
0(+  t)  = y < 1 and  o ' (^  t)  =0; 

and 

(2.4)  is  nonlattice  . 

Define  the  random  function  X by 

n 
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«*•«**  rntinn.  . 


X (t) 
n / 


!]nt] 

an* 


ostsi, 


r* 

where  of  = Q"(x)/y,  0 < a < ®,  and  [x]  is  the  greatest  integer  in  x. 


Our  goal  in  this  chapter  ii  to  prove  that  the  finite-dimensional 

distributions  (f.d.d.’s)  of  the  random  function  X , conditioned  on 

n 

n < N < “,  cor; /erge  weakly  to  those  of  Brownian  excursion,  W*. 

To  be  more  specific,  we  assume  that  k i 1}  are  the  coordinate 

functions  defined  on  the  product  space 


(n,  j,  p)  = X (R> 


n=l 

where  R = (-«>,  <») , 31  is  the  o-field  of  Borel  sets  of  R,  and  n is 

the  common  probability  measure  of  the  £ *s.  If  A = {r  < N < °°),  then 

we  let  (A  > A A P ) be  the  trace  of  (H,  S,  P)  on  A , where  A H J 
n n n ' 7 ' n n 

= {An  fl  F,  F e f)  and  Pn(A)  = P(a)  pfA^)  for  A e Afl  fl  SF. 


Our  result  has  an  application  to  queueing  theory,  as  well  as  the 

obvious  interpretation  as  the  fortune  of  a gambler  or  insurance  company 

prior  to  ruin.  If  is  the  waiting  time  of  nth  customer  in  a general 

single  server  queue,  then  the  £'s  in  this  application  are  differences 

of  service  and  interarrival  times.  In  this  context  N is  the  number 

of  customers  served  in  the  first  busy  period.  Observe  that  W = S 

n n 

on  {N  > n}.  Thus  conditioning  on  {n  < N < »}  will  yield  limit  theorems 
for  the  waiting  time  process,  given  that  the  first  busy  period  has  not 
ended  but  will  end  eventually. 

We  close  this  section  by  noting  a well  known  fact  that  P(N  < »}  = 1 
when  < 0 [see  CHUNG  (1968)  p.  244],  In  this  case  we  can  write 
{N  > n)  instead  of  {n  < N < »}. 


2.  Some  Asymptotic  Results 

Conditions  (2.1)  through  (2.4)  imply  that  as  n -* 


(2.5) 


>{Sn  £ 0)  ~ (2jtn)“^  A»0"\  if  ux  < 0, 


and 


(2.6)  P{Sn  ^ 0}  ~ ( 2rtn)  2 Y (cyt)  , if  ^ > 0. 


These  results  are  contained  in  BAHADUR  and  RAO  (1960);  also  see  IGLEHART 
(1974b). 

Our  analysis  relies  heavily  on  two  lemmas,  which  appear  in  IGLEHART 
(1974b). 


(2.7)  LEMMA  [IGLEHART  (1974b),  Lemma  2.l],  Let 


w | w 

Z d/  = exP  ) z vn 


n=o 


n=l 


for  |/6|  S 1.  If  b = 0(n  2), 


then 


d = 0(n  ' } as  n -*  «>. 


(2.8)  LEMMA  [IGLEHART  (1974b),  Lemma  2.2].  Let  c , d SO,  c ~ cn~^, 

E d n d < ® and  d = 0(n-1).  If  a = c d , then 
n=o  n n — n j=o  n-j  y 


a /N'  cdn 
n 


•4 


as  n -»  co. 


Let  fn  = P{N  = n)  and  r = P{n  < N < ®).  Our  first  result  yields 
the  asymptotic  behavior  of  the  sequence  of  [f^:  n s 1}  and  (r^  n s 1). 
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(2.9)  THEOREM.  If.  conditions  (2.1)  - (2.4)  hold,  then  as  n -+  » 


(2.10) 


V1 


A 


n ~ T / « 


n2  (2*)V 


(2.11) 


7 A 

cn  ~ “3 "X — ~IT > 

n n2  (2Jt)5ctT(7  1 - 1) 


where 


A = 


exp 


09  y~n 

- E — p(s  * 

' n 11  n 


n=l 


(Y  1 -l)  exp 


E xl 

n=l 


V[Sn  > 0} 


if 


> °; 


if  p^  < 0. 


Proof . When  p^  < 0,  the  asymptotic  relation  (2.11)  is  proved  in 

IGLEHART  (1974b).  Theorem  2.1.  Because  f = r - r , it  follows 

n n— l n 

immediately  that  given  g > 0,  there  exists  nQ  such  that  for  all 


n > n. 


f ? 


Y 


+ °(®)  * fn  s -^—4 

Y -1  n2 


+ o(e) , 


where  f = (2rt)  ^(<yx)  ^A.  Since  g is  arbitrary,  we  obtain  (2.10)  for 

U-f  < 0. 

To  obtain  (2,10)  and  (2,11)  for  p^  > 0,  we  use  cn  identity  from 
random  walk  theory  [see  CHUNG  (1968),  p.  256-258], 

1 - E ^*n  f exp(itS  )dP  = exp [-  E ~ f expfitS  )dp] 

(N=n)  " [ £l  n (8^0}  n'  J 
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for  m < 1.  Setting  t = 0 yields 


(2.12) 


-E 

n=l 


f h - exp 
n 


2 — P{S  £ 0) 


n=l 


n n 


Differentiating  both  sides  of  (2.12)  with  respect  to  A yields 


(2. 


00 

13)  V nf  T"n>in  = exp 

5=1  n 


oo  ..-n 

E — pts„s0)i” 

n=l 


oo 

Y'  Y'r‘P{Sn£0}/% 


n=l 


where  | < 1.  Set  an  = n^v"  , cn  = 7~  P{Sn  £ 0] 


and 


E V"  * *4 E V '«•. * °^”1 

n=o  L n=l  J 


Then  from  (2.13),  we  see  that  an  = *j|£  Cn-jdj*  N°W  Cn  ~^2,tn)  (aT)  ’ 


because  of  (2.6).  So  b^  = n cn  = 0(n  2).  Also, 


E - 


= exp 


n=o 


00  y~  n 

-E  V PiSn  S °> 

n=l 


< 00 . 


Hence  we  can  apply  Lemmas  2.7  and  2.8  to  obtain  (2.10). 

Now  for  any  e > 0,  there  exists  nQ  such  that  for  all  n a nQ, 


(2.14) 


yfl  y1 

(i  - e)f  -T  £ f_  * (1  + e)f  — 

n2  n2 


Hence 
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£ (i  + e)f 


n 


3 

2 


Vn 

Y r„ 


£ 


(l  + e)f 


CO 


k=n+l 


yn+k 


00 


E t 

k=l 


k 


Letting  n -»  ®,  we  obtain 


(2.15)  lim  sup  n1  7~n±  £ (l  + SyfCT”1  -l)"1  . 

n-Ko  n 

Similarly,  using  the  first  inequality  of  (2.14),  we  have 

(2.16)  lim  inf  n%  T~nf  ^ (l  - e)^"1  -l)"1  . 

n->  <»  n 

Combining  (2.15)  and  (2.16),  we  obtain  (2.11),  since  e is  arbitrary. 

Our  problem  is  closely  related  to  the  asymptotic  analysis  of  the 

distribution  function  of  M , where  M = maxfS,  : 0 ^ k £ n).  This 

subject  has  been  treated  extensively  by  BOROVKOV  (1965).  The  following 

two  theorems  are  immediate  results  of  BOROVKOV  (1965),  Theorem  7. 

Denote  M = sup{S  : k S 0). 

K 

(2.17)  THEOREM  . .If.  conditions  (2.1)  - (2. 4)  hold,  then  for  x > G as 
n -*  o° 
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r*atrojK^ 


0 ; 


« i 


m 


ls&; 

& 


wv*WVW  •'A  M 


(2.19) 


( Mn  ) l H 

P{ £ x > -P 


- £ xj  ~ ~~  exp(-xa Vo)  A 1 } 

n ) tC  cn 


lu/n  ' (o/n 

The  limits  in  (2.18)  and  (2.19)  are  uniform  for  x in  compact  sets 


for  (j,  < 0. 


(2.20)  THEOREM.  If  conditions  (2.i)  - (2.4)  hold . then  for  x > 0 as 


n -»  « 


(2.21) 


and 


(2.22) 


Pin  < N /—  < 00 


-xaTn  ^ J 


T T:v— 

■n  ot-x  { y - i, 


for  (A,  > 0, 


pKwr>  ni  ~ r TtTr"  ■ 


QfT  (Y  - 1) 


for  [A.  < 0, 


The  limits  in  (2.21)  and  (2.22)  are  uniform  for  x _in  compact  sets . 

2 , 

Recall  that  R(x)  = x expf-x  /2),  x £ 0,  is  the  density  of  the 
Raleigh  distribution.  Theorem  2.17  enables  us  to  understand  the  asymp- 
totic behavior  of  S , conditioned  on  fN  > nl  as  n -*  «>. 

(2.23)  THEOREM  . If  conditions  (2.1)  - (2.4)  hold . then  for  x > 0 as 
n -»  os 

| — — £ x,  N > n|  ~ cxp(x&V»)  > if  |i1  > 0, 


<2-24>  I aJS 


n at 


and 
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(2.25) 


> 


if  |i  < 0 


N > n 


fw 


y1  r(x) 

— expC-xaVn)  ~1 — 

it2  at 


The  limits  in  (2.24)  and  (2.25)  are  uniform  for  x in  compact  sets. 

Proof.  We  shall  prove  (2,24)  first.  Recall  that  = Wn  on  {N  > n). 
So  if  we  can  show  (2.24)  with  Sn  replaced  by  W^,  the  result  will  be 
established.  A simple  path  decomposition  yields 


' W 


(2.26) 


\afn 


£ x,  N > n 


1 = p 

Vn  ^ 

| - E p| 

} 

(a/n 

k=l  ( 

W 


£ p{~  5 x>  N = k|  . 

o/n 


Using  the  fact  that  N is  an  optional  random  variable  and  that  WN  = 0, 
we  can  rewrite  (2.2v)  as 


(W 

(2.27)  P -2-  <;  x,  N > n 
iotjn 


(W  ) 

P — s x 
(o/n  ) 


n 


E 

k=l 


V 


(w  ^ 
| n-k 

io/n 


For  n ^ 1 and  x S 0,  let 

-n  (W 

h(n,x)  = n7  exp(  -xon+/n)  P\ £ x 

^<*/n 

For  e > 0 and  a fixed  0 < n^  < [n(l  - c)]  we  can  write,  because  of 


(2.27), 
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(2.29) 


lira  h(n,x) 

ir*» 


Mk1s-1 

jr2crr 


and 


(2.30) 


lira  I 
n-x»  n 


R(x}  a-i 
j^xr 


as- 


_ ,;2/2 

because  of  (2.18).  Since  R(x)  = xe  '■  is  bounded,  relation  (2.18) 

also  implies  that  h(n,x)  s>  G (0  < G < »)  for  all  n s 1 and  x <:  0.  Using 
this  fact,  we  can  conclude  chat 

tn(l-e)] 


«■»’  - * x * • *-1  j; 


-k 

T fk  • 


k=n  +1 
o 


Because  of  (2.10),  y ”fn  s Hn  2 (0  < H < ®)  for  nil.  Thus  we  have 
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HP***1- 'ytS» m •■— -■—  | '^..,-y.-^~.^ 


K £ GII 
n 


k=[n(l-e)  ]-f-J. 


;:  GH 


E (.-»)• 


[n(l  - €)]2  k=[n( 1-e) J+l 


Finally, 


e)"2  n 


_3  1 

£ GH(l  - e)  2 n *(l  + In  e + In  n)  . 


(2.33) 


, . -n  -xaWn  . _ 

lim  nT  e v f =0 

]»CO  n 


because  of  (2.10)  and  being  positive.  Using  (2.28)  - (2.33),  and 
selecting  n^  sufficiently  large  we  can  make 


«r" 


arbitrarily  small  for  large  n.  Since 


A = exp 


[-E^ks^o]].  ,-E 


-k  . 

y f. 


by  virtue  of  (2.12),  this  completes  the  proof  of  (2,24).  Because  (2.24) 
is  a direct  consequence  of  (2.18),  that  the  limit  in  (2.24)  is  uniform  for 
x in  compact  sets  follows  immediately  from  Theorem  2.17. 
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To  establish  (2.25)  for  the  case  < 0,  we  can  write  from  (2,26) 


ILl 

(o/n 


) 

W ) 

n 

w , 1 

x,  N>n > = P 

— > x> 

- £ f.p< 

n-k  ^ ( 

> x> 

) 

r 

Ot/n  j 

6i  k 

o/n  j 

{M_  ) 

( M ) 

pj-2-  > x j 

- Py ^ x/ 

(o/n  J 

(o/n  ) 

n 

Ev 

k=l 


M , ) 1 M 

P<_£±  > x ( - P<— — > x! 


<*/n  j (o/n 


+ pi—  > x 1 P{N  > nj 
{a/n  ‘ 


M 


p{-a-*  J - pU-*  xi 


~r  l " T7 

o/n  J (o/n 


n 


+ E «„ 

k=]  K 


Jjsi  SJ  Jji-J 

|a/!T  J («£  j 


+ pj—  >x|P{N>n} 
(o/n  ) 


For  x £ 0 and  n 5 1,  let 


g(n,x)  = n7  e}:p(xQfX>/n) 


JX  £ 4 -pU;  £ x! 

) Wn  J 


!a/n 


Then  from  (2.19) 


(2.35) 


lim  g(n,x) 

n->co 


R(x)  A-l_ 


Multiplying  both  sides  of  (2.34)  by  n7  n exp(xo;x Vn) , we  have 


(2.36) 


nT  n exp(xaT\/n)  p|-^cr  > x,  N > n> 

|a/n  ' ( 

n-l  . — r— 

= ~g(n, x)  + V r"  f • g(n-k,  xjn/n-k)  . 

k=l  k 

- n7  " exp(xaVn)  P|j^  > x j p{n  > n ). 


To  estimate  the  last  terra  appearing  in  (2.36),  let  K be  the  nontrivial 
solution  of  0(/6)  = 1 . Then  for  x > 0 as  n -»  oo 


(2.37) 


P{  M > xo/n)  ~ D exp(  -xotnjn) 


where  D is  a positive  constant  whose  precise  value  need  not  concern 
us;  see  IGLEHART  (1972),  Lemma  1 for  this  result.  Hence, 

(2.38)  lim  n7  ° exp(xoctjn)  p|-^->x(  p{N  > n)  = 0 

n_>”  l«/n  ) 


because  of  (2.11),  (2.37),  and  K > t >0.  The  summation  term  in  (2.36) 
can  be  broken  into  three  parts  and  be  taken  care  of  in  the  same  manner  as 
we  have  done  in  (2.28)  — (2.32).  And  we  can  conclude  this  time  that 


(2.39) 


Umron7^xp(xc*x^)  p|^=  > x,  N > nj 


co  n 

A 


The  power  series  ^ P{Sn  > 0)  has  > 1 as  radius  of  converg- 
es? — n 

ence  and  E PfS  > 0}  is  finite;  see  VERAVERBEKE  and  TEUGELS 

n=l  n ' n 1 

(1975),  p.  281,  for  a discussion  of  this  result.  Therefore,  from  (2.12) 


we  have 


( 2.40) 


-i  + E T'”i  » (f'-o  exprf£! 


e*P  £ — P{s„  > 0] 


n n 


Combining  (2.39)  and  (2.40)  we  obtain  (2,25).  Uniform  convergence  is  a 
direct  consequence  of  Theorem  2.17.  This  completes  our  proof. 


For  x £ 0 and  n £ 1,  we  define 


(2.41)  F 


n(x)  - (nit)^  >'  j p)~~;2dy,  N>  n|  exp(+yotv/n)  , 
(o,xJ  (cWsAi  ) 


where  the  sign  is  taken  for  u > 0,  and  the  V sigi)  is  taken 


for  p1  < 0. 


(2.42)  corollary,  j^r  ail  x ;>  o, 


lim  F (x)  = 1 

n-K»  n 


the  Raleigh  distribution. 


<4), 


_Proof.  Integration  by  parts  yields 


Fn(x)  - (njr)^7  p|-^:  S x,  N > n | exp(-x«x/n) 


+ n at  n7 


'*  f p)^- 

(o,xj  lotjn 


s y.  N > n>exp(-yor&yn)dy  , 


if  u1  > 0.  Hence,  by  Theorem  2.23,  we  have 


Fn(x)  = °(D  - ] y exp/-.  |-)dy 
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If  jj^  < °,  integration  by  parts  will  yield 


Fn(x) 


i _T  <w  I 

njt)  T P{n  > n}  - Pj— 2—  > x N > n>  exp(xat/n) 

<c*/£  ) 

+ reform  )'  n ^ pj_IL_  > y}  N > n | exp(ya'W^i)d 
(o,xj  la/n  l 


N > n>  exp(yartv/n)dy 


= o(!)  , \ y exp(-  ^-) 


(o,xj 


because  of  (2.11)  and  Theorem  2.23  . 


The  next  . suit  is  needed  in  proving  the  convergence  of  high  dimen- 
sional distributions. 


(2.43)  THEOREM.  Assume  conditions  (2.1)  - (2.4)  are  satisfied.  For 


x,  y > 0 as  n -» 


(2.44) 


— y 

oan2 


exp{(x  ■ yjort/n}  g(l,  y,  x). 


where 


g(t>y.x)  = (2ift)  11  |c\p  j - exp|^-  j. 


The  convergence  in  (2.44)  uniform  for  x compact  sets . 


Proof.  Let  M = max{-S  : 0 £ k ^ n}  and  M = sup{-S  :k  a 0).  If 


we  set 


then 


Ifn 
(a/n 


a„(y,  x)  = pj-il-  £ x - y,  N_ya/~  > 


(M  (-S  ) 

an(y,  *)  = p H:<  y, 


(c*/n 


a/n 


a y - x>  . 


Using  the  simple  identity  P(A  fi  B)  = P(A)  + P(A  Pi  B)  - P(B),  we  can  write 


(M"  ) (Mn  (-S)  ) {(~S  ) 

.(y,*)  - PH=  < y + p{~  » y>  — r~<  y~*i  -p{— ■£-<  y-x 


(2.45) 


\ctfi  yi  l<*/5  ^ <*/£ 


ot/n 


= I + J - K 
n n n 


But  we  have 


(2.46) 


xn  = k yj  + °(x)  exp(-y(>Wn)  , 


[<</£ 


(2.47)  J 


(2.48)  K = 
n 


= |p(“L.<  y] ^r-  exp[-(y-x)crrN/n-^~^-l|[l  + °(l)] 

n | I (2«j»)*aT  L 2 J) 

ji £-  exp T- (y-x) - ■ ( [l+ o(l)]  ; 

( (2«n)*ar  L 2 


see  BOROVKOV  (1965),  Theorem  7,  for  the  above  results.  Combining  (2.45) 
to  (2.48),  we  obtain  (2.4-1).  Uniform  convergence  also  follow*  directly 
from  BOROVKOV  (1965;,  Theorem  7. 


3.  The  Main  Limit  Theorems 

Having  obtained  various  asymptotic  results  in  the  last  section,  we 

are  now  able  to  establish  the  main  result  of  this  chapter.  We  shall  show 

the  convergence  of  one-dimensional  distribution  of  (X  n < N <.  ®),  then 

n ' 

extend  the  result  to  higher  dimensions.  First,  we  state  a standard  result 
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in  the  theory  of  weak  convergence  used  by  IGLEHART  (1974a),  Lemma  (2.18). 
Alternatively,  see  BELKIN  (1972),  p.  54,  or  BILLINGSLEY  (1968),  Theorem 
5.5  . 

(2.49)  LEMMA.  Let  {p^:  n S 1}  be  a sequence  of  finite  measures  on  ft, 
the  Borel  sets  of  R = (-»,  + ®) . Suppose  p^  =?>  p finite.  If  {f^: 
n £ 1 ) _i_s  a_  sequence  of  uniformly  bounded . Borel  measurable  functions 
converging  uniformly  on  compact  sets  to  an  everywhere  bounded  continuous 
limit  f , then 

lim  ) f (x)p  (dx) 
n^oo  j'  n n 

for  B e 9.  provided  p(bB)  = 0. 

The  following  two  results  are  the  main  objective  of  this  chapter. 


(2.50) 

THEOREM. 

If  conditions  (2.1)  - 

(2.4)  are 

satisfied,  then  for  all 

x £ 0 

and 

0 £ t 

SI  as  n -»  °> 

(2.51) 

( X (t)  | n < N < «) 

=$>  wj(t) 

• 

Proof. 

The 

claim 

(2.51)  is  trivial  for 

t = 0. 

Next,  consider  t = 1, 

and  take  0 

< a < 

b.  We  want  to  estimate 

( S 

(2.52)  P a < 5 b | n < N < 03 

l <Vn 


\ f(x)p(dx) 
B 


which  is  equal  to 


(2.53) 


r 1 P l a < —2—  s b . n < N < » ! 

n ' Q^r  \ 


- r_1  J p|_JL_  edx  N > nj 
n (a,b]  1«A  1 


pxQ(/n  ^ < oo  j . 


If  p > 0,  then  (2.52)  is  furthermore  equal  to 


(2.54) 


( itn' 


“T — j (nit)^7"n-  exp( -xar/n) ■ pf-^-  6dx,  N > nl 

)*rn  (a)b]  lo/IT  ’ 

• n exp(x<Wn)  P[N  x(^r  < °°} 


■f  C 

— — T — I F (dx)  • n exp(xd?tv/n)  pfN  r < <»}  . 

-M-  (A  i n'  1 “xcn/n  J 


(*n  )5r  (a,  b] 


Now 


(2.55)  n exp(xofrfn)  p(N_xQ^  < °°}  - n exp(xon\/n)  P{M  > xos/n} 

~ D • n expf-xaC^-TVn  ]. 

Combining  (2.11),  (2.52)  - (2.55),  and  using  Corollary  2.42  and  Lemma 
2.49,  we  conclude  that  for  > 0 


(2.56) 


•USl. 


lim  P a < -2-  £ b n < N < 00  f = 0. 

l 0^ 


If 


H < 0,  then  (2.52)  is  less  than  or  equal  to 


-1  \ ° 

r P<  - ~ ■ •>  a.  N > n 

n Wn 


which  is  asymptotically  as  n -♦  ® equal  to 


^2^  exp( -xOftJn)  R(x)  • (r  ^ -l)A 


because  of  (2.11.)  and  (2.25).  This  shows  that  (2.56)  also  holds  when 


< 0. 


-&5- 


- -v^f-  2B°«*s£2tSS«*9ff^?3£^^A3 


Next,  consider  0 < t < 1 and  take  0 < a < b.  For  all  choices  of 
a,  b,  and  t it  will  suffice  to  show  that 


(2.57) 

If  n1  > 0, 

(2.58) 

is  equal  to 


lim  P{a  < 
n-»°° 


lint] 

<Vn 


£ b n < N < » 


?{ a < W*(t)  £ b}  . 


a < JnSII  ^ b j n < N < oo 

*JZ 


\ , f pi^  ®dx,  N > [nt]jpXQs/“  £n-[nt]  < N < ~} 
( a,  b]  la/ri  ) 

- - “j  ! j (nt*)fy  ntexp(-xor^i) 
[«t(l-t)n3]srn  (a/b] 


.(flsil,  * 


l<K/St  Ji 


e — • , N > [nt] 


n(l-t)r  1 t^exp(xo'Vn)pjn-[nt ] < N . 


Appealing  to  Corollary  2.42,  we  see  that  (2.57)  equals  further  to 

(2.59)  -f-’—T- f Ffnt](T>-^ir|^>(*W^)P{n-[nt]<lI_^<»). 

n rn  T w 

Because  of  (2.21), 


(2.60) 


A 


Combining  (2.11),  (2.41),  (2.59),  (2,60)  and  using  Lemma  2.49,  we  obtain 


A 


(**>  (1'*>  (a,bl^ 


<w(  y"  -1)  (T-l)  Jt 


2x  exp(~ 


2 1 


rbtj)dx 


37 


[2ntJ(i-t) 


1 


(a,b] 


Pq(0,0,1:,x)  dx  ( 


which  proves  (2.57)  for  p > 0. 

Similarly,  if  p < 0,  the  term  (2.58)  becomes 

P ja  < -tSil  <:  b | N > n]  , 

( cm/ n ) 

which  is  equal  to 


(2.61) 


y1 


i 

(rtt)^(i-t) 


C /dx\ / _ ( ) 

) F[nt](^V(X^)  exP(-x^)p|N.XQ^>n-[nt)j 
(a,b]  V 


Using  the  same  argument,  we  see (2. 61)  converge  as  n -»  » to 


j P*(o,  0,  t,  x)dx  . 
(a,  b] 

This  completes  the  proof. 


Finally,  we  can  plunge  into  the  proof  of  convergence  of  the  f.d.d.’s 


(2.62)  THEOREM.  If  conditions  (2.1)  - (2.4)  hold,  then  for  k > 1 

and  0 s t ,<  t.  < ...  <t,  5 1 as  n *♦  ® • 

12  k — 

(2.63)  ...,  Xn(tR)  | n<N<«)  =*>  (wj(t,),  ...,  wj(tk))  . 


From  (2,60)  we  see  that  the  quotient  term  in  the  last  product  term  is 
asymptotically  equal  to 


(2.65) 


\i-t  / 

' m ' 


y exp 
m 

2 *1 

ym 

y . exp 

r 2 

- ym-l 

m-i 

2(l- 1 

u Hh* 

expr(Wi)£W^i  • 


Let 


fin  (S[nt  l-[nt  3 

G-<  y~ ) = 3 p{ cdx,  N ..  r> rntl-tnt  ] 

1 -y  -ttyo  m m-i. 

m-1 


n m ' 


0+  ( Oi/n 


-y£n*nJ  Lntm_^]  exp|^  (x-y^  ^ | . 
Integration  by  parts  yields 


fm 

<2-66)  - 5 «(*.  - w Vi-  x)dx' 

because  of  Theorem  (2,43).  Finally,  by  our  induction  assumption 


(2.67) 


lim  P 


[ntl] 

Oifn 


£ X^j  • • •> 


’[nt  l3 
m-1 

Ctyn 


^ x n < N < a# , 
m-1  1 


x x , 
•1  em-1 


. r •$ 


0 J0 

dym-l*'*  dyi  • 


Combining  (2  . 64)  - (2.67)  and  using  Lemma  (.1,49)  twice,  we  obtain 


29' 
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CHAPTER  3 

CONDITIONED  ' IMIT  THEOREMS  FOR  RANDOM  WALKS 
WITH  POSITIVE  DRIFT 


1.  Introduction 

Let  (X  : k i 1)  be  a sequence  of  independent  random  variables,  iden- 
k 

tic-.-.ily  distributed  with  common  distribution  F(x)  on  (-oo,  +oo).  Let 

2 2 2 

EjX^)  = and  EfX^  - = a , 0 < cf  < °o.  Form  the  random  walk 

fS  : n S 0)  by  setting  S = 0 and  S = X + •••  + X , n S 1.  Next 
1 n 0 nl  n 

let  N be  the  hitting  time  of  the  set  (-oo,  x]  by  the  random  walk 
x 

N = inf  fn  > 0:  S ^ x}  . 
x n 

where  the  infimum  of  the  empty  set  is  taken  to  be  +oo.  Set  N = Nq. 

When  p^  > 0,  the  random  walk  is  drifting  to  +oo  and  asymptotically 
the  conditioning  N > n plays  no  role:  the  random  walk  does  not  feel 

the  barrier  at  the  origin.  In  fact,  we  have 

lim  p|-^-  S x | N > nj  = 0,  x>0 

IW“  W } 

because  of  Theorem  2.23  and  P{N  = w)  > 0 [see  CHUNG  (1968),  proof  of 
Theorem  8.4.4],  In  this  case  it  is  more  meaningful  to  examine,  instead, 
the  limit  distribution  of  S^,  conditioned  on  (n  < N < oo).  From  Theorem 
2.50  we  have,  for  0 £ t £ 1, 

lim  pj-Iaii  £ X I n < N < ®)  = P{w+(t)  £ x]  . 

< o/n  > 
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But  P{W^(1)  <:  x)  = 1 and  this  result  does  not  tell  us  the  limiting  be- 
havior of  Sn,  because  normalization  by  n^  kills  it.  It  suggests  that 

1 

we  need  not  normalize  Sn  by  nB  to  get  convergence  to  a non-degenerate 
limit. 

The  main  result  in  Section  2 is  to  show  (S  n < N < oo)  converges 

n 

A 

to  a non-degenerate  random  variable  S , and  to  identify  the  distribution 
A 

function  of  S , As  we  have  pointed  out  before,  this  problem  is  closely 
related  to  the  asymptotic  analysis  of  the  distribution  function  of  Mr, 
where  = max{S^:  0 £ k £ n).  We  shall  get  the  asymptotic  result  of 
P{Mn  £ x)  in  Section  3. 

2.  Limit  Theorem  For  (S  ^ | n < N < oo) 

Throughout  this  chapter,  we  shall  assume  that  the  distribution  of 
X1  satisfies  the  following  conditions: 

(3.1)  0 < £ oo; 


(3.2)  0(sj  = E exp(sX^)  converges  for  real  s g (-a,0],  for  some 
a > 0; 


(3.3)  6(s)  attains  its  infimum  at  a point  -t,  0 < r < a,  where 

6(t)  = y < 1,  and  0*(t)  = 0; 


and 

(3.4)  if  Xx  is  lattice,  then  Pj^  = 0}  > 0 . 

The  proof  of  the  main  result  in  this  section  requires  us  to  intro- 
duce the  so-called  associated  random  variable  to  X^;  see  FELLER  (1971), 
p.  406,  for  a discussion  of  this  concept. 


The  distribution  function  F(x)  of  the  associated  random  variable 
is  given  by 

rx 

t(x)  = | exp(-Ky)  F(dy)  . 

iCO 

Let  {X  , 5L,  ...}  be  a sequence  of  i.i.d.  copies  of  X . Form  the 
2 3 

associated  random  walk  {£  : n s 0}  by  setting  S - 0 and  S = 

n u n 

$+...+X,n£l.  Also,  define  the  stopping  time  N accordingly.  If 
§(.4)  = E expC&ft.)  is  the  moment  generating  function  of  X^,  then 
§C6)  = 6 (A  - K).  Note  that  0'(-k)  < 0 implies  that  the  associated 
random  walk  {§  : n £ 0}  has  a negative  drift.  Furthermore,  the  distribu- 
tion  function  £ of  5^  satisfies  conditions  (1.7)  - (1.10).  Hence,  we 
have  from  (1.14) 


(3.5) 


lim  p{§  ^ x j N > n}  = c • v(x)  , 


where  c is  a constant  and  £(x)  is  a known  function  solely  depending 

on  x;  see  Ch.  1,  Sec.  2,  for  this  result. 

To  get  the  limit  result  for  (Sn  | n < N < co),  we  first  estimate 

PfS  £ x,  n < N < oo}  for  x > 0. 
v n 

(3.6)  PfS  ^ x,  n < N < «}  = j P[S  6 dy,  N > n}py{N  < oo} 

n / J n n 

(o;xj 


= I(o  x](yi  + "*+  yn^dyl^"'F^dyn^  P^M  > yl  + "*+yJ 


where  a =(( y, , . . yr ) : y.>  0,  i£  k < n}  . Exploiting  the  concept 


of  associated  random  variable  to  X^,  we  can  write  (3.6)  as 

r"|l  (y,  + ••*+y„)  p[m“  >y  +**‘+y  }exp[K(y  +-..+y  )Jp(dy  )...p(dy  ) 
J M in  1 n 1 n 1 n 


[ P{M"  > y]eKy  P{&  £ dy,  N > nj  . 

(o;x] 


Hence 


P(s  £ x,  n < N < °o} 
P{S  £ x|n  < N < <»}  = 2 

n P{n  < N < oo} 


(o, x ]P^M~  > y}e*yp£sn  S dy|N  > n} 
f p{m‘  > y}e<ypfs  6 dy|N  > n} 

(o;oo]  * 

It  follows  from  (3.6)  and  Lemma  2.49  that 

( P[U~>y}eKy  V(dy) 

(3.7)  lim  P{S  ^ x|n  < N <S»]  = ■ H(x)  . 

11-+®  n f 

r Ai  P{M“>y}e^(dy) 


It  is  easy  to  see  that  H(x)  is  nondecreasing  and  right  continuous  with 
H(0)  = 0,  H(-h»)  = 1.  Thus  H(x)  is  a distribution  function. 

Hence,  we  arrive  at  our  main  result  in  this  section: 


(3,8)  THEOREM.  Jjf  conditions  (3.1)  - (3.4)  are  satisfied,  then  as  n -»  oo 


(3.9) 


(S  | n < N < »)  =$>  , 


1 s a non-degenerate  random  variable. 
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One  application  of  this  result  is  in  terms  of  the  waiting  time  of  the 

nth  customer,  W , in  a general  single  server  queue  with  traffic  intensity 

p > 1.  In  this  case,  W = S on  the  set  {N  > n)  and  N is  the  number 

n n 

£ 

of  customers  served  in  the  first  busy  period.  Thus  S can  be  thought 
of  as  the  limiting  waiting  time  given  that  the  first  busy  period  has  not 
ended  yet  but  will  eventually  terminate. 

For  illustration,  we  take  as  example  X,  = v-  u,  where  v(u)  has 
an  exponential  distribution  with  parameter  p(\).  This  corresponds  to 
M/M/1  queue.  If  P = l/n,  we  need  to  assume  P>  1 in  order  to  insure 
Ex^  = p^  > 0.  Then,  it  is  well-known  that 

(3.10)  P{M  > y)  = p 1 exp(-  Ky)  ; 

see,  for  example,  FELLER  (1971),  p,  199.  Because  of  (3.5),  (3.7), 

(3.10)  and  V(oo)  = c 1,  we  have 

H(x)  = cV(x)  , 


which  is  the  distribution  function  of  lim(S  ft  > n).  One  c .1  calculate 

n-»oo  n ’ 

its  Laplace  transform  and  find 

(3.11)  E exp(uS^)  = J-L-.R}-iu.+A]  , 

X(2u  + X - ,■) 

see  IOLEHART  (1974b),  p.  750,  for  details. 

3.  The  Asymptotic  Analysis  of  the  Distribution  Function  of  M 

Following  the  exposition  of  the  last  section,  we  shall  go  on  to  estab- 
lish the  following 


-35- 


(3.12)  THEOREM.  If  conditions  (3.1)  - (3.4)  are  satisfied,  then  for  all 
x :>  0 


(3.13) 


p(Mn  * 


x) 


V1  . G(x 


2 

n2 


(2it)2qt 


exp 


-k 


Z -V  p(sk  * °} 

k=l  k k 


as  n oo,  where  G(x)  is  a known  function  solely  depending  on  x. 

Proof.  For  x £ 0 let  G (x)  = PfM  £ x]  and  u (x)  = PfS  £ x.  N > nt. 
n n n 1 n 

n & 1.  Employing  the  same  argument  that  is  used  in  deriving  (3.7),  we 
obtain 

u (x)  = j eKy  P(S  e dy,  N > n) 

(o,xj 

= P{Sn  5 x»  N > n}eKX  - k J P(Sn  S y,  b n}e^ydy. 


The  last  equality  is  due  to  integration  by  parts.  VERAVERBEKE  and  TEUGEES 
(1975)  have  shown  that 


as  n -*  oo.  Hence 


Pfs  s£  x,  N > n) 
v n ’ 1 


/ V(x) 
n2  (2it)^orT 


(3.14) 


l,n(x) 


/ V(x) 
n1  (2r)2aT 


where  V\x)  = V(x)e*X 
for  {un(x):  n S 0 } : 


K,  J V(y)ekydy. 
(o,  x J 


Form  the  generating  function 


= £ u (x)/jn  . 

n=l 
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Then  we  have  the  following  Spitze.  Jentity 


(3.15)  Gn(x)-4n  = [l  + U(A,x)  J • exp 

n=l 


oo  n I 

E irp(s„ s °) 

n=l 


(see  [30],  P.  280).  Differentiating  with  respect  to  A yields 


£n7~nG  (x)/ 
n=l 


(3.16) 


= ~ir  p^sn  5 0)/j\ 

•fjl]  7'nP{SnS0}/j[x+u(|,x)]  + Sn7"nun(x)/ 

I (n=l  ' ' n=l 


1 


As  we  had  in  (2.6), 


(3.17)  P(Sn  * 0}  - (2-n)"^7n(aT)"  as  n -»  00 . 


If  we  set 


00 

“ 

(3.18) 

a (x)  = 

n 

Z T'nP{sn  * 0}.:n 

1 

\l(~,  x) 

[_n=l  j 

> 

- 

then,  because  of  (3.14),  (3.17)  and  Lemma  2.8, 


a (x)  ~ a(x;n  . as  n -*  cc  } 

n 


where  a(x) 


(SsJ'^aT)'1  [1  + U(7_1,  x)  j.  Thus 


(3.19) 


a (x)  i n7 


fx) 


~ (2r>n) 


U(7  1,x)  r V(x) ] 


TT>  .J.V  l ,1. 


by  virtue  of  (3,14).  Let 


(3.20) 


o°  r co  -n  *1 

E V"  - exPl  T P£Sn  * °)^ 

n=o  Ln=l  J 


Then,  because  of  (3.17)  and  Lemma  2.7, 


(3.21) 


dn  = &[n~z)> 


as  n *»  oo  , 


And  it  is  obvious  that 


KS?pSiJfce?S6SSi’!?5iS  aiKS= 


(3.22)  d < co  . 

uj  n 
n=o 

Combining  (3.16),  (3.18),  and  (3.20),  we  have 


n-o 


n7~nG  (x)/ 
n 


00 

T d 

n 

• < 

* 

E [ 

_n=o 

l n=o  1 

a (x)  + nY  ”u  (x)l/5n 
n n J 


Because  of  (3.19)  - (3.22)  and  Lemma  2.8,  we  obtain 


nY  ~ ( 2nn)  ^(cr.)  1G(x)  • exp 


“ y~k 

S -TT  p<sk  s °) 


k=l 


where  G(x)  = 1 + U(Y  \ x)  + V(x).  This  completes  our  proof. 
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CHAPTER  4 

CONDITIONED  FUNCTIONAL  CENTRAL  LIMIT  THEOREMS  FOR  QUEUES 
WITH  TRAFFIC  INTENSITY  EQUAL  TO  UNITY 


1.  Introduction 

We  shall  first  study  a general  single-server  queueing  system  with 
traffic  intensity,  p,  equal  to  1,  and  then  excend  the  result  to 
multiple-channel  queues.  The  Gl/G/l  queueing  system  is  constructed 
in  Section  1 of  Ch.  1.  We  shall  follow  the  exposition  there. 

Section  2 of  this  chapter  is  concerned  with  studying  the  behavior 
of  the  workload  process,  W(t),  conditioned  on  the  fact  that  the  first 
busy  period  has  not  ended  by  time  t,  as  t becomes  large.  A func- 
tional central  limit  theorem  in  D[0,  1]  is  our  goal.  In  Section  3, 
the  queue-length  process,  Q(t),  is  treated. 

In  Section  4 we  shall  extend  the  results  for  the  W and  W(t) 

n 

processes  to  a variation  of  Gl/G/l  queueing  system,  where  there  are 
several  initial  customers  instead  of  only  one  at  time  t = 0.  This  simple 
generalization  enables  us  to  establish  results  for  our  multiple-channel 
queueing  systems,  which  is  the  subject  of  Section  5. 

2,  The  Workload  Process  fW(t):  t a 0) 

For  t S 0 let  A(t)  = n+1  or.  ft  £ t < t 1,  where  as  previously 

n n+l 

defined  in  Section  1 of  Ch.  1,  t = u^  + +un»  n a 1,  and  tQ  = 0. 

For  convenience  let  A(O-)  = 0.  Clearly(  fA(t):  t a 0}  is  a renewal 
process  which  represents  the  number  of  arrivals  in  the  interval  [0,  t]. 
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Set  I(t)  = vQ  + •••  + i - t.  The  workload  process,  W(t),  is 


represented  as 


(4.1) 


W(t)  = I(t)  - inf {l( 4) : 0<  & < t}  . 


with  W(0+)  = vQ  > 0,  Define  the  random  function 


v*)  ■ » 


os  t <1 


Next,  define  a stopping  time  L = inf{t  > 0:  W(t)  = 0),  where  L repre- 
sents the  length  of  the  first  busy  period.  Our  goal  in  this  section  iB 
to  obc«.in  a central  limit  theorem  for  (Yy|L  > u)  as  u -*  ».  The  first 
step  toward  this  goal  is  to  understand  the  asymptotic  behavior  of 
du  = P{L  > u}. 


(4,2)  LEMMA,  For  u > 0,  P{L  > u}  ~ c(\  u)  ® j|S  u -►  »,  where  c is 
a constant. 


Proof : 


(4.3) 


{L  > u)  = (W(t)  > 0,  0 < t < u} 


= {W(t1)  > 0,W(t2)  > 0, ...,W(A(u)-l)  > 0,  W(u)>0) 
- <Sl>0-—SA(u)-l>  °-  SA(U)+  *A(u)-U>0)  ’ 
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and 

(4.4)  {SjL  > 0:  1 < i < A(u)J  C {L  < u)  C {Si  > 0:  1<  i < A(u)-l) 
It  has  been  shown  that 

(4.5)  P{S1  >0:  1 < i < A(u)3  ~ c(\u)~^,  as  u -»  co; 

see  IGLEHART  (1974a),  Lemma  4.2,  The  desired  result  follows  from  (4.4) 
and  (4.5). 

The  proof  of  the  main  theorem  requires  a standard  result  in  the 
theory  of  weak  convergence.  Let  {X^:  n £ 1}  and  {Y^:  n £ 1 } be  two 
sequences  of  random  elements  of  a separable  metrix  space  such  that  X 

n 

and  have  a common  domain.  Let  p(x,y)  = sup{  |x(,6)  - yCA)|  :0  S ^ S 1} 
for  any  x,  y e D[0,  1], 

(4.6)  THEOREM.  li_  Xfi  =$>X  arid  p(*n,  Y ) =$>0,  then  Yr  =$>X. 

See  BILLINGSLEY  (1968),  Theorem  4.1,  for  this  result. 

Let 

zU)  = , 0 < t < 1 . 

o(Xu)^ 

2 

(4.7)  LEMMA.  _If  p = 1,  and  0 < 0 < to,  then  as  u -»  <» 

(Zj  L > u)  =5>  W+  . 
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Proof . It  follows  from  Lemma  4.2  and  IGLEHART  (1974a),  Lemma  4.3|that, 
if  0 < t,  < t„  < ...  < t S 1, 

(^y(  tj)  ) • • •)  Zu(  tj^)  |L  > u)  =J>  (W  (t^))  . 

Thus  we  have  the  convergence  of  f.d.d.'s  of  (Zu|l  > u).  It  remains  to 
show  that  the  family  ( (zu | L > u):  u £ 0}  is  tight.  Since  zu(°)  = 0> 
it  suffices  to  show  that,  for  every  6 > 0, 

(4.8)  lim  11m  p[co  (8)  > e | L > u}  =0, 

610  U400  zu 

where  aWS)  = sup  { |x(s)-x(t)|:  s,  t e [0,  1],  |s-t|  < 6};  see 
BILLINGSLEY  (1968),  Theorem  15.5,  for  this  result. 

We  have  from  (4.4) 

0 £ P(u>  (8)  > e I L > u) 
z 

u 

< £fN  > Mu).:.1)  (8)  > e | N > A(u)  - 1}  . 

P{L  > u)  u 

1GI£HART  (1974a),  Lemma  4.6,  has  shown  that 

lim  lini  P{o>  (8)  > e | N > A(u)  - l)  =0. 

6lo  u^oo  z„ 


Hence  (4.8)  follows  immediately.  This  completes  the  proof. 
Now  we  are  able  to  establish  our  main  result. 


(4.9)  THEOREM,  If  p = 1,  0 < cr  <oo  and  E{  u ^ ) < oo(  then  as  u -* 


(Y„  | L > u)  W+ 


Proof.  On  the  set  (L  > u), 


A(u)-1 

W(u)  = £ v - u 


SA(u)  + (tA(u)  ' 


Let  p(x,y)  = sup{|x(t)  - y(t)j:  OSti  1)  for  any  x,  y in  D[0,  l]. 


P<Y„'  SUPC 1 ^(u^)  " UZI=  0i/^l) 


cj(Xu) 


t sup{u  , s : 0 < A < l} 

cr(Xu)® 


T sup{ u ; 1 < k < A(u)}  . 

cr(\u)* 


Thus  P{p(Y  , Z ) > e L > u}  is  less  than  or  equal  to 


i? 


-1 

d 

u 


p;  sup 
tl<k<A(u)+l 


u. 


> e cr(Xu) 


) 


< 


;a(u)+i 

s 


1{uk>ea(\u)^ ) 


d"1  E 

U 


|a(u)+i|  • pju^  > g ct(Au)^^ 


The  last  equality  holds  because  of  Wald’s  identity.  It  follows  from 
Lemma  4,2  that 


(4.10) 


p(p(v  zu)  > « I L > u) 


= [l  + o(u)  ] 


-1  E{A(u)+1) 

a 


u 


u2P{u1>ecr(A.u)*} . 


The  elementary  renewal  theorem  says  that 


(4.11) 


lim 

U-*00 


E(A(u)3 

u 


\ . 


Our  assumption  E(u^)  < oo  implies  that 


(4.12) 


lim  t3P{u  > t)  =0  . 

t-*>o  1 
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Combining  (4.10)  - (4,12),  wa  have 


(4.13) 


u»p{p(V  z.)  > e I l > u|  =0 


Finally,  the  result  follows  immediately  from  (4.6),  (4.7),  and  (4.13). 

i. 

Now  let  us  suppose  that  at  time  t there  are  rQ(t)  - 1]  cus- 
tomers waiting  in  the  queue  and  the  server  is  serving  some  customer,  say 
j*"*1  customer,  whose  residual  service  time  v*  is  part  of  v^ . Then, 
we  can  rewrite  the  workload  process,  W(t),  in  the  following  way: 


W(t)  = 


A(t) 

E 

i=A( t)-[Q( t)-l ]+ 


v.  + v' 
x 


i.e.,  the  sum  of  [Q(t)  - 1]+  complete  service  times  and  a residual 
service  time. 

Le^  V(t)  = W(t)  - v'.  Since  the  difference  of  V(t)  and  W(t) 
is  dominated  by  max{v  : 0 £ k £ A(t)},  it  is  easy  to  see  that  V(t), 

K 

properly  normalized,  has  the  same  weak  limit  (see  the  proof  of  Theorem 
4.9).  For  0 £ t s 1,  let  V^(t)  = V(ut)/c(xu)^. 

(4.14)  THEOREM.  _If  p = 1,  0 < cf2  < » and  Efu^+v^3)  < ®,  then  as 

u -»  oo 

(vu|l  > u)  =>  w + . 
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3,  The  Queue-Length  Process  fQ(t):  t a 0] 

We  now  turn  to  the  queue-length  process,  (Q(t):  t & 0],  which 
represents  the  number  of  customers  in  the  system  at  time  t,  including 
the  one  being  served.  This  process  is  a random  step  function  on  the 
positive  half  line  with  jumps  of  + 1.  We  may  define  Q(t)  formally 
by  introducing  the  departure  epochs  {d^}: 


[t  + v + W , if  n > 1, 
l n n n'  ’ 


if  n =0. 


The  departure  process  {D(t):  t & 0),  which  records  the  number  of  depar- 
ture in  (0,  t],  is  defined  as 


D(t)  = 


max(n>l  : d St),  if  d < t 
n ’ 1 ’ 


, if  dx  > t. 


Clearly  Q(t)  = A(t)  - D(t),  where  of  course  A(t)  and  D(t)  are  highly 
dependent.  We  observe  that  A(t)  ^ D(t)  for  all  t s 0 and 


(4.15) 


(N>  D(  t) ) Z>  {N>  A(  t) ) , for  all  t > 0. 


Let  us  consider  another  renewal  process  {D'(t>:  t 0). 


D’(t) 


max{n£l:  v +»»*+vn£t),  if  v^  < t, 


, if  v L > t. 


Note  ,’r.at  D’(t)  and  they  are  equal  if  there  are  no  idle  periods 
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in  (0,t],  We  observe  that  the  renewal  process  (D’(t):  t ^ 0}  has 
the  same  rate  as  {A(t):  t ^ 0}  because  p = 1.  Using  the  same  argument 
as  in  IGLEHART  (1973),  Lemma  4.2,  we  have 

(4.16)  LEMMA.  P(N>D*(u)}~  c(m;u)-^  as  u -*  “,  where  c is  _a  constant . 
The  precise  value  of  c is  not  our  concern  here;  see  FELLER  (1971) 

p.  415,  for  its  value. 

Jk 

(4.17)  COROLLARY.  P{N  > D(u) } ~ c(p-u)  as  u ®. 

Using  relation  (4.15),  (4.17)  and  the  same  prooi  employed  in 
IGIEHART  (1973),  Lemma  (4.2),  we  have  the  following  two  results. 

2 3 3 

(4.18)  LEMMA.  If.  p = 1,  0 < c <co  and  E{u  + v ) < then  _as 

t -»  CD 

/ | N > D(t)\  =*>  W+(l)  , 

\ff(*t)*  / 

(4.19)  COROLLARY.  If  p = 1,  0 < o2  < »,  and  Efu^+Vj3}  < °o,  then  for 
x > 0 

lim  pjv(t)  < x | N > D( t) | = 0. 


Corollary  4.19  enables  us  to  show  that,  as  time  goes  to  infinity,  the 
queue  length  will  also  become  large,  conditioned  on  the  event  that  the 
first  busy  period  never  ends.  For  the  unconditioned  case,  a similar 
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i 

: \ 

? 

i 

I 


result  is  first  proved  in  WHITT  (1968),  pp.  149-151. 

(4.20)  LEMMA.  For  any  M > 0, 

lim  P{Q( t)  < M | N > A( t) } = 0 . 

t->co 


Proof.  Let  (A^,  At  H 5,  p^)  be  the  trace  of  (ft,  5,  P)  on  A = (N  > D(t)}, 
We  can  write,  for  any  x and  M, 


P {V(t)  < x}  = £ P (V(t)  < x,  [Q( t)  - lj 4 = k}  . 
k=o 


Let 


E = { V(  t)  < x}  = 


_ I 


A(t) 


I 


( i^(tj^Q(t)-ll+  1 j' 


ms 


F = C Cq( t ) - if  = k}  , 

and  let  sj(x^,  ...,  x^)  denote  the  sigma-field  generated  by  random  vari- 
ables x^,  x^.  Then  on  the  set  (A(t)  = i+k,  D(l)  = i},  for  any 

i,  k > 0, 

E e ’<Vi’  •••>  Vk5  ' 


^ ® °(V  Vj/  Ui>  "*>  °i+k^  ’ 


and 


A t e 3(vQ,  vi_i,  uL,  ...,  uA)  . 
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Therefore,  p(Ef|A  ) = P(e)  • P^F^), 


i.e. , 


Pt(V(t)  < x)  = 


A(t) 

s £ 

k=o  |i=[A(t)-k]+ 


Ti 5 xr p 


- k} 


M ( A(t)  ) 

> p < vi  - 4-  = k| 

k=0  (i=[A(t)-k]+  j 


For  any  M,  choose  x sufficiently  large  so  that  for  all  k ^ M, 


( A(t) 

! 1 

A(t)  ■ 

£ + ^x 

f S ^ 

£ 

|i=[A(t)-k]  ^ 

i-[A(t)-M'J 

The  choice  of  x is  independent  of  t,  because  the  last  inequality 
involves  only  at  most  m v^s,  which  are  independent  and  identically 
distributed.  Thus  we  have 

Pt(V(t)  < x}  > e Pt([Q(t)  -if  < M}  . 

Since  the  left  side  goes  to  zero  as  t goes  to  infinity  by  Corollary 
4.19,  so  does  the  right  side.  The  theorem  follows  from  the  inequality 


p{n>d( t) } 

P(Q(t)  < M|N  > A(t))  < P.(Q(t)  < M). 

P{N>A(t)) 


Using  essentially  the  same  proof  employed  in  (4.20)  we  have  the 
following  generalization. 
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(4.21)  COROLLARY.  For  an£  M > 0 and  0 < t * 1, 

lim  P{Q( ut)  < M j N > A(u)}  = 0. 

Vt*co 


Our  main  object  of  interest  is 


Q„(t) 


Q(at)  . 

7(?uj5  ' 


0 < t < 1 . 


We  are  interested  in  the  limiting  behavior  of  (Q^Jl  u).  The  next 
lemma  is  helpful  in  establishing  the  convergence  of  the  f.d.d.’s  of 
(QjL  > u). 


(4.22)  LEMMA.  For  any  e > 0 and  0 < t < 1, 


y£p|®$'x"11  26 1 L>uj  = °- 

Proof . Set  E(u)  = (k  : |—  - X | < g)  and  F(u)  its  complement.  Let 
Pu(a)  = P(a|n  > A(u))  for  A e (N  > A( u) } n J and  r = P{N  > A(u)).  We 
shall  first  consider 

<4-23>  &]■*  @ -1'1' 2e , * >*<“>!• 


Probability  (4.23)  is  less  than  or  equal  to 


- \ 1|  > e,  A(ut)€E(ut),  D(ut)eE(ut),  Q(ut)  > M| 

+ Pu[Q(ut)  < M}  + r“1|p{A(ut)€F(ut)}+P(D,(ut)eF(ut)}j 


As  u •»  “,  Ku  is  asymptotically  equal  to 


£uip(A(ut)  c F(ut»  < ) _ 

6 t U2 


which  will  converge  to  zero  because  the  numerator  is  finite;  see  SMITH 
(1958),  pp.  248-249.  Because  of  Corollary  4.21,  also  converges  to 
zero.  It  remains  to  estimate  1^.  Because  of  (4.15), 


I < r"1  . pj|  ~ 

u u , I' 

k,feE(ut)  [ 


E (v.-X'1) 
.=£+2 

k - l 


Since 


A(ut)  = k,  D( ut)  = l,  N > & j . 

jA(ut)  = k,  D(ut)  = l,  N > i!  j eo(v0»  v£+1>  •••>  uk+1) 


( E (vi 
/ I i=i+2 


k - £ 


* - j €0<W 


'**>  vk)  , 


the  two  events  are  independent.  It  follows 


I i 
u 


r . p; > € .• 

u i-  if.  tr/ui ) I k - l 1 


k,  ikE1 

k-£>M 


(ut)  | 


p{k>£,  A(ut)=k,  D(ut)=£} 


By  Chebyshev’s  inequality  we  have 
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i , ^51^2  £ 


u u 


2 

€ • M 


k,  ikE(ut) 
k-i>M 


p[n  >[n(X-e)  ],  A(ut)=k>  D(ut)=0} 


e{(v1-\"1)2|  p|n  > [n(X-e)]| 


Putting  the  estimates  of  Iu>  Ju  and  together,  and  letting  u go 
to  infinity  in  (4.23),  we  have 


f|V(ut)  ) E{(v1-X"1)2}  h 

lim  sup  P <1  •■-■'  ■ - X ^ e?  ~ 2 * * ( \ ) • 

iwoc  u(Q(ut)  i €jm  * 


Since  M is  arbitrary,  we  have 


!V(ut)  , 

|— ; r - X > e |N  > A(u) 

Q(utj 


Using  the  basic  relation  (4,4),  we  have  the  result. 


(4.24)  COROLLARY.  For  any  £ > 0 and_  0 < ^ < • • • < tfe  s 1 , 


( v(utJ  i ) 

-”1W'1"l>£;l£1£k|L>ul 


Proof.  Since 


P< 


V(ut . ) 
i 


( Q(ut.) 


- X I > e;  1 < i < k L > u 


k ( V(ut  ) j 

5 L P)H — r - x i > e | l > u , 

i^i  ( Q(ut1)  ) 


result  follows  immediately  from  Lemma  4.22. 

Now,  we  are  in  the  position  to  show  the  convergence  of  the  f.d.d.'s 
of  (qJl  > u). 

(4.25)  THEOREM.  If  p =■■  1,  0 < a2  < co 
0 < t.  <•«■•<  t,  ^ 1,  then  as  u ~*  oo 

(QyC^),  ...,  QU(tk)|L  > “) 


and  Efu^+v^3}  < co,  and  if 

=>  (w+(t1);...,w+(tk))< 


Proof.  First,  suppose  k = 1.  We  shall  show  that  for  any  x > 0 and 
0 < t £ 1, 


(4.26) 


(Q(ut)  ) , 

lim  P< ^ — t-  < x I L > u}  = p{w  ( t)  < x} 

»***  |cr(X3u)4  ) 


Let  P^(A)=  p(A|L  > u)  for  any  A e {L  > u}  (1  J,  Then 
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|Q(ul)  ) _ (v(ut)  v(ut)  ^ ^ ) 

? < T < X>  = P \ T £ Xx> 

u(tf(X3u)2  ) u(cr(Xu)2  Q(ut)  ) 


s p )v(ut^  , v(utl  t Xx>  _ x_1|  < e 


♦ ',S  - 5 «! 


< p jllylL.  < (i  + xe)xl  + 


e, 


because  of  Lemma  4.22.  Since  e is  arbitrary 


1 Q(uti_  < xl  < 


(4.27)  lim  sup  P t < x 

u->»  u|  \i 

ValX  u)2 


i 


lim  p 
u-»» 


W * 


= P{W+(t)<x} 


On  the  other  hand. 


( Q(ut)  1 

P“WSX< 


( Q(ut) 

> p < — - — r - x, 

uia(X3u)* 


V(ut) 

Q(ut) 


X | > e 


( v(ut)  ) 

2:  p < r < (l  - eX)x>  . 


u(a(Xu) 


T 


Again,  because  e i f»  arbitrary, 
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(4.28) 


lim  inf  P 
u *»  » 


| Q(ut) 


< x;  > lim  P 


u-*» 


u 


V(ut) 


cr(xu)w 


S X, 


Equality  in  (4.26)  follows  from  (4.27)  and  (4.28), 

Exactly  the  same  argument  can  he  employed  to  establish  the  result 
for  k > 1.  This  time  we  use  Corollary  4.24  instead.  And  Theorem  4.25 
follows. 

In  order  to  get  the  full  weak  convergence  in  E>[0,  1],  we  need 
to  show  that  the  family  {(Qy|L  > u):  u ^ 0}  is  tight.  To  do  so  it 
suffices  to  show  that  for  any  e > 0, 


(4.29) 


lim  lim  P<o^  (5,  0,  l)  > € I N > A(u)>  = 0 . 

5io  u-»®  ( Qu  ' ’ 


see  BILLINGSLEY  (1968),  Theorems  15.1  and  15.5  for  this  result. 

The  first  step  in  demonstrating  (4.29)  is  the  next  lemma. 


(4.30)  LEMMA.  For  every  e > 0, 


(4.31) 


lim 

Tlo 


lim  r 
vwcc  u 


p 


sup  ■ 

,o£A£r  a(X  u)2 


> 


e. 


N > A(ut) 


0 . 


On  the  set  {N  > A(ut)),  (Q(it6):  0 ^ A ^ t)  is  a stationary 

process  with  independent  increment.  In  fact,  Q(u,6)  = A(uA)  - D*(a4), 

0 £ h £ t,  i.e.,  the  difference  of  two  standard  renewal  processes.  Based 
on  this  observation,  the  proof  of  Lemma  4.30  follows  the  argument  of 
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BELKIN  (1972),  p.  49,  and  will  be  omitted.  Also  see  IGLEHART  (1973), 
pp.  19-24. 

With  this  lemma  in  hand  it  is  an  easy  matter  to  show 


(4.32)  THEOREM.  If  p = 1,  0 < cj2  < co  and  Efu^+v^}  < oo,  then 
(Q  Jl  > u)  =5>W+,  as  u -»  oo. 


Proof.  As  remarked  above,  it  suffices  to  show  (4.29).  For  every  x e(0,l], 
® > 0 and  0 < 6 < x 


(4.33)  pCw  (», 0,1)26}  < r"plo>  (8,0, l)  > e,  N > a(ux) 

I yu 


S r“lp|\(8'0>l)iS'  ^ < «.  «(») 


Q(a-i) 


+ r P{  sup 


sup  — r 

U ( o<^><t  ij(X3u)^ 


> €,  N > A(ux) 


= I + J 
u u 


Now.  I is  bounded  above  by 
' u 


(4.34) 


r”1  P]<h  (8,  -r-8,  l)  > e,  N > A(ux)j 
u ' ^ u 


£ r,’,1  PK  (%  1)  - N > D(u( t-8) ) 


On  the  set  {A(u(x-6))=  k,  D(u(x-6))=  i) , k > i ;>  0,  we  see  that  two 


events 
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{N  > d(u(t-8)  )}  e tf(vQ,  ...jV^;  ...,  u^) 


{m  (8,T-8,l)  > e}  e cr(v,  , •••>  uk+1>  •••) 
wu 

are  independent.  Thus,  (4.34)  is  further  less  than  or  equal  to 


P{o>Q  (8,t-S,i)  > e] 
Qu 


P{N  > d(u(t-8))} 


P[%  (5,0,1'/  > e} 
wu 


Since  Qy  =£>  f(B)  and  P{f(B)  e c}  = 1 where  B is  the  standard 
Brownian  motion  in  [0,  l],  it  follows 


(4.35) 


lim  lim  PjV  (8,  0,  l)  > e}  =0 
ft-K>  U->“  W,, 


see  WHITT  (1968),  Theorem  8.11,  for  this  result.  Combining  (4.31)  - (4.35), 
we  obtain  (4.29).  This  completes  the  proof. 


4.  A Single  Server  System  With  Several  Initial  Customers 

In  this  section  we  shall  consider  a simple  variation  of  the  standard 

Gl/G/l  queue  with  p = 1.  We  assume  that  at  time  0,  instead  of  one 

initial  customer,  there  are  m (m  > 1)  initial  customers.  Thus  our 

queueing  system  is  defined  by  the  two  sequences  [u^:  n 2:  1 ) and  [v^: 

n a - m + 1},  where  u^  represents  the  interarrival  times  between  the 

nth  and  (n+l)St  customers,  and  v the  service  time  of  nth  customer. 

n 
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Let  X = v - u , n S 1.  Form  the  random  walk  fS  : n > 0}  by 
n n-m  n’  1 n ' 

setting  S = 0 and  S = X,  + • • . + X f n s 1.  And  define  the  other 
0 n 1 n 

quantities,  e.g.  N,  L,  W^,  W(t),  accordingly.  Then  the  next 

theorem  is  an  immediate  result  of  Theorems  1.6  and  4.9. 


(4.36)  THEOREM.  t£  p = 1, 


0 < a <oo  and  E{  u J } < 


then 


(4.37) 


as  n -»  co 


We  observe  that  N represents  the  total  number  of  customers 
being  served  during  the  first  period  for  which  there  are  always  at  least 
m customers  in  the  system.  The  purpose  of  deriving  (4.37)  and  (4.38) 
is  that  these  results  will  be  helpful  in  establishing  conditioned  weak 
limit  theorems  in  multiple  channel  queueing  systems.  This  is  dealt  with 
in  the  next  section. 


5.  Multiple  Channel  Queues  With  Identical  Servers 

In  this  section  conditioned  limit  theorems  for  queues  with  traffic 
intensity  equal  unity  are  extended  to  multiple  channel  queueing  systems. 
The  systems  are  defined  by  two  basic  sequences  of  independent  and  identi- 
cally distributed  random  variables  {u^:  n s 1]  and  (v^:;  n ^ - m + 1), 
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atilt  aitteiMifc «■»*»>•■ 


where  represents  the  interarrival  time  between  the  and  (n+l)st 

customers,  and  v the  service  duration  of  nth  customer.  There  are  m 
identical  service  channels  and  one  initial  customer  for  each  service 
channel.  The  queueing  discipline  is  first  come  first  served. 

Let  Efu^)  = A 1 and  E{v^)  = p,  \ We  assume  that  p = A/ny,  = 1. 
Let  A(t)  be  the  arrival  process  associated  with  {u^:  n s 1).  We  shall 
define  the  gross  input  G(t)  and  net  input  I(t)  of  the  systems  as 
follows:  for  t ^ 0 

a(  t) 

G(t)  = S Vn  ’ 

n=-m+l 

I(t)  = G(  t ) - mt  . 

Let  H(t)  denote  the  total  workload  of  the  systems  at 
Define  the  workload  at  time  t in  the  jth  channel,  H^(t), 
total  work  that  is  eventually  done  by  server  j in  time  t. 

H( t)  = H, ( t)  + ...  + H (t),  t > 0. 

1 m 

If  we  define  the  virtual  waiting  time  W(t)  as 

w(t)  = min{H.(t):  1 r j < m] , t > 0, 

3 

then  W(t)  represents  the  time  that  a hypothetical  customer  arriving  at 
time  t would  have  'o  .ait  before  being  served  by  the  first  available 
server.  Define  a stopping  time 

L = inf { t > 0:  w(t)  = 0)  . 


t ime  t. . 

to  be  the 
So 
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So  before  time  L all  the  servers  are  busy. 

For  the  moment  consider  a single  server  system  A with  m initial 
customers.  System  A is  described  by  the  basic  information  {(W  , u^)  = 
n & 1],  where  u’  = u , v'  = v /m  and  the  traffic  intensity  is  unity. 
From  a workload  point  of  view,  the  multiple  channel  systems  behave  exactly 
like  system  A on  the  set  {L  > u).  Based  on  this  observation,  we  have 
the  following  result  by  virtue  of  (4.38): 


(4.39) 


where 


(jM 

\c(Xu)* 


L > L =5>  W , 


c = 


’ 2,  , a <vl> 

cr  (u2)  + 


m -> 


On  the  set  {L  > u}  for  any  i,  .]  and  0 i.  t s;  1,  H^(ut)  and 
Hj(ut)  may  differ  at  most  by  the  service  duration  of  one  customer  who 
arrives  before  time  u,  that  is  to  say,  for  0 <.  t <,  1, 


|H.(ut)  - H.(ut)|  < max{v  : - m+1  < k < A(t)} 

i j k 

Following  the  argument  used  in  proving  (4.13),  we  see  that  for  every 
e > o, 


(4.40) 


1 

( /Hi<u,) 

H^(u.) 

\ , 

lim  P 

IWOO 

(PVC(Xu)^ 

c(Xu)^ 

• > e 1 l > u 

= 0 . 
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iwiim  a 1 1 


\ip 


:kv'j 
v.«  ■ ,*-- 


,;V, 


Hence  all  processes 


'»!<“•>  | \ 

| L>  uL  1 £ i < m , 


c(Xu) 


¥ 


have  the  same  weak  limit,  if  there  is  one.  But  H(t) 
Using  the  triangular  inequality,  we  have 


H^t)  + 


(4.41) 


H(u*  ) x ^ /Hk^U*)  H-:(U*) 


~X  > 


cm(Xu)2  c(Xu)^/  m k^j  \c(Xu)5  c(Xu) 


\¥*  777~\¥ 


Combining  (4.39)  - (4.41),  we  obtain 


(4.42) 


Hj(-) 


cm  *(Xu) 


\ + 

— r | L > u =5»W,  as  u -)  m 

.u)4  / 


+ H (t). 
m 


for  any  j = 1,  2,  .. . , m. 

Since  W(t)  = min[H^(t):  1 £ j ^ m},  it  is  easy  to  see  that 


(4.43) 


/ m 

p(w(u.),  H^uOj  < Y,  p(VU’)'  Hj(u*)) 

k=l 


Combining  (4.42)  and  (4.43)  and  appealing  to  Theorem  4.6,  we  establish 
the  following. 


X 2 3 

(4.44)  THEOREM.  If  0 = ~ = 1,  0 < a < « and  Ef  ^ )<  °°  then 


-2^L,  |L>u\ 

cm  (Xu)*  j 


W 


as  u -*  oo  , 
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•ft*  #*fcvw  *t*^>*m**XM^^tt**L. % 


where 


2 «2(v.)li 

- ”<V+V-  • 

%—  — 


We  shall  complete  this  section  with  a weak  limit  theorem  for  the 


process  (W  : nil},  the  waiting  time  of  n^h  customer.  Note  that 


can  be  obtained  by  the  random  scaling 


(4.45) 


w = w(t  ) ; 

n n ; 


where  t = u + •••  +u  . Define 
n 1 n 


N = sup{n  > 0:  w(4)  > 0 for  all  A £ t } 


where  the  supremum  of  the  empty  set  is  taken  to  be  + m.  Then  we  have 


the  following 


2 3 

(4.46)  THEOREM . _If  p = 1,  0 < 0 <oo  and  E{  u ) < oo,  then  as 


n -»  oo, 


Wr  , \ 

--Is-U  , \ + 

-1  I N > n / =$>  w . 
cm  n / 


Proof . First  we  define  a random  change  of  time  e^  by 


5^(4)  = min^-^^-  , 1^  , 0 < 4 £ 


"K^»  - W[a6]  . 


But  W(^)  = S'  + (t  - A)  on  the  set  E , 
A(A)  AU)  n’ 

Sn  =kh  (vk  " uk}*  Hence 


where 


W 


(4.47) 


-1 

cm  n‘ 


A(nenU)) 

7TJ — 


cm  n 


+ e U) 

n 


> 


on  the  set  E , where 
n 

(4.48)  (e  j N > n)  =>  0 . 

n 

We  can  assume  Efu^}  = X 1 < 1 (if  this  were  not  the  case,  a 

__  n 

linear  rescaling  of  time  will  make  \ x < 1).  Using  Chebyshev's  inequality, 
we  can  easily  show  that  for  0 £ A £ 1 


(4.49) 


lim  P1 
n-»“ 


Ikil>  x | 


N > n 


= 0 . 


Combining  (4.37),  (4.47)  to  (4.49),  and  Theorem  4.6,  we  prove  the  result. 
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sequence  of  independent,  identically  distributed  random  variables.  Form 

the  random  walk  fS  ••  n s 0]  by  setting  = 0 and  S = t + ...  + t 
1 n 0 n^l  cn 

n fe  1,  and  let  N denote  the  hitting  time  of  the  set  (~»,  o]  by  the 
random  walk:  N = inf{n  >0:  s 0}.  In  queueing  context  N represents  j 

the  number  of  customers  served  in  the  first  busy  period.  The  natural  j 

measure  of  congestion  for  the  system  is  the  traffic  intensity  p = \/\i, 
where  \ and  p,  are  the  arrival  and  service  rates,  respectively. 

The  basic  processes  we  are  interested  in  are  {W^:  n £ 0},  ! 

{W(t):  t £ 0},  and  {Q(t):  t ;>  0}  where  represents  the  waiting  time  j 

of  nth  customer,  W(t)  denotes  the  workload  of  the  server  at  time  t,  and  j 
Q(t)  is  the  number  of  customers  in  the  system  at  time  t.  Define  the  ; 

random  function  Xr  by  <t ) = where  a is  a norming  constant,  j 
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When  p /£  1 and  appropriate  conditions  on  the  distribution  of  ^ j 

hold,  it  is  shown  that  the  finite-dimensional  distributions  of  X , con- 
ditioned  on  n < N < »,  converge  weakly  to  those  of  Brownian  excursion, 
as  n goes  to  infinity.  When  p > 1,  it  is  shown  that  W , conditioned 

on  n < N < »,  converges  weakly  to  a limit  random  variable  S*,  and  the  j 

M 1 

distribution  function  of  S is  exhibited.  ■ 

j 

These  problems  are  closely  related  to  the  asymptotic  analysis  of  the  I 
distribution  of  the  maximum  of  the  random  walk  {S^:  n 2 o).  Various  j 

asymptotic  results  on  the  joint  distribution  function  cf  the  maximum  and  !■ 

N are  obtained,  ! 

If  p = 1 and  some  mild  regularity  conditions  hold,  the  workload 
process,  W(t),  suitably  normed  and  conditioned  on  the  event  that  the 
first  busy  period  has  not  terminated  yet,  is  shown  to  converge  weakly 
to  Brownian  meander,  as  time  goes  to  infinity.  The  same  result  for  the 
queue-length  process,  Q(t),  is  also  obtained.  These  results  are  extended 
to  a variation  of  the  general  single  server  queueing  system,  where  there 
are  several  initial  customers  instead  of  only  one.  This  simple  generaliza- 
tion enaoles  us  to  establish  results  for  certain  multiple-channel  queueing 
systems. 
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